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ABSTRACT 

In this paper we develop a consistent, phenomenological methodology to measure 
preferred-frame effects (PFEs) in binary pulsars that exhibit a high rate of periastron 
advance. We show that in these systems the existence of a preferred frame for gravity 
leads to an observable characteristic 'signature' in the timing data, which uniquely 
identifies this effect. We expand the standard Damour-Deruelle timing formula to 
incorporate this 'signature' and show how this new PFE timing model can be used to 
either measure or constrain the parameters related to a violation of the local Lorentz 
invariance of gravity in the strong internal fields of neutron stars. In particular, we 
demonstrate that in the presence of PFEs we expect a set of the new timing parameters 
to have a unique relationship that can be measured and tested incontrovertibly. This 
new methodology is applied to the Double Pulsar, which turns out to be the ideal 
test system for this kind of experiments. The currently available dataset allows us 
only to study the impact of PFEs on the orbital precession rate, ti;, providing limits 
that are, at the moment, clearly less stringent than existing limits on PFE strong- 
field parameters. However, simulations show that the constraints improve fast in the 
coming years, allowing us to study all new PFE timing parameters and to check for 
the unique relationship between them. Finally, we show how a combination of several 
suitable systems in a PFE antenna array, expected to be availabe for instance with 
the Square-Kilometre- Array (SKA), provides full sensitivity to possible violations of 
local Lorentz invariance in strong gravitational fields in all directions of the sky. This 
PFE antenna array may eventually allow us to determine the direction of a preferred 
frame should it exist. 
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1 INTRODUCTION 

The theory of general relativity (GR) has so far passed all ex- 
perimental tests with flying colours (Stairs 2003, Will 2006, 
Kramer et al. 2006). Nevertheless, GR may not be the flnal 
word in our understanding of gravity and indeed, alternative 
theories of gravity exist, predicting deviations from GR in 
various possible ways. 

Some theories of gravity predict that the Universe's 
global matter distribution selects a preferred rest frame for 
local gravitational physics. In such theories the outcome of 
gravitational experiments depends on the motion of the lab- 
oratory with respect to this preferred frame. In particular, 
theories in which gravity is partially mediated by a vec- 
tor field or a second tensor field are known to exhibit such 
preferred-frame effects whose strength is determined by cos- 
mological matching parameters. 
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In the post-Newtonian limit, such preferred-frame ef- 
fects are described by two phenomenological parameters, the 
parameterised post-Newtonian (PPN) parameters oi and 
02. In GR these two parameters are zero. Solar system ex- 
periments have already tightly constrained these two PPN 
parameters. For details see Will (1993, 2006) and references 
therein. On the other hand, limits obtained in the weak grav- 
itational fields of the solar system cannot rule out effects 
that only become significant in strong gravitational fields. 
For the well-motivated scalar-tensor theories, this has been 
demonstrated in a series of papers by Damour and Esposito- 
Farese (1992a, 1993, 1996a, 1996b). 

In order to describe preferred-frame effects in binary 
pulsars, Damour and Esposito-Farese (1992b) introduced 
a non-boost-invariant Lagrangian that includes two strong 
field parameters 6t\ and 0.2 , capable of accounting for strong- 
field violations of a local Lorentz invariance of gravity. In the 
weak-field limit these two strong-field parameter are equal to 
Qi and Qf2, respectively. Concentrating on the contribution 
by Qi, Damour and Esposito-Farese show that qi has the ef- 
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feet of inducing secular variations in the orbital eccentricity 
and in the orientation of the pulsar orbit, i.e. the location 
of periastron and the orientation of the orbital plane. These 
variations depend on the magnitude and direction of the 
velocity with respect to the preferred frame. Their calcula- 
tions for small-eccentricity orbits shows that the eccentricity 
oscillates in a specific way between a maximum and a min- 
imum value, with the period of the relativistic precession 
of periastron. FYom this, using probabilistic considerations, 
they wore able to derive an upper limit of |cii| < 5 X 10"'' 
with a 90% confidence level. After the discovery of a num- 
ber of new small-eccentricity binaries, Wex (2000) extended 
this method in statistically combining multiple systems, by 
this taking care of a potential selection effect when simply 
picking the system with the most favourable parameter com- 
bination. His analysis yielded the slightly improved limit of 
|di| < 1.2 X 10""* with a 95% confidence level. 

In a recent paper Bailey and Kostelecky (2006) study 
the violation of local Lorentz invariance within an ofl^ective 
field theory called the standard model extension (SME). Be- 
sides a number of solar system tests, they use their formalism 
to investigate the measurability of SME preferred-frame ef- 
fects in binary pulsars. However, as already pointed out by 
Bailey and Kostelecky, their formalism does not account for 
strong field effects. In addition, as we will demonstrate in 
the course of this paper, their usage of the PSR B1913-I-16 
timing results does not provide a consistent test. 

In this paper we present a new method for testing 
preferred-frame effects related to strong gravitational fields, 
which uses binary pulsars with a high rate of periastron 
advance. Our method provides independent and simultane- 
ous tests for both di and a2. While the tests of Damour 
and Esposito-Farese (1992b) and Wex (2000), by their prob- 
abilistic considerations, can only provide upper limits for 
preferred-frame effects, the method introduced in this pa- 
per is even capable of detecting a preferred frame in the 
universe, if it exists. 

We apply our new method to the PSR .J0737-3039 bi- 
nary pulsar system (Burgay et al. 2003; Lyne et al. 2004), 
which is not only the only Double Pulsar system known, 
but which also exhibits the highest rate of periastron ad- 
vance, LO = dio/dt, of any binary pulsar known. The 
precisely measured value of w = 16.89947(68) deg yr~^ 
(Kramer et al. 2006) is more than four times larger than the 
value measured for the Hulse- Taylor pulsar PSR B1913-t-16 
(Weisbcrg & Taylor 2002). While the Double Pulsar has 
been used to provide the so far most stringent test of GR 
in the strong field regime, we show here that the large pe- 
riastron advance makes it also a unique tcstbed to actually 
measure a possible violation of the Lorentz invariance for 
strong gravitational fields, which may occur in alternative 
theories of gravity. 

The plan of the paper is as follows. In Section 2, based 
on a theory-independent framework of Will (1993), we study 
in full the motion of a compact binary system in theories of 
gravity that are non-boost-invariant. We calculate the sec- 
ular evolution, caused by preferred-frame effects, for all bi- 
nary parameters, in order to show that the combination of 
these periodic parameter changes exhibits a unique 'signa- 
ture' of a preferred frame. In Section 3 we investigate the im- 
pact of such time-dependencies on the high-precision timing 
measurements of binary pulsars, and how this can be used 



to either measure or constrain strong-field preferred-frame 
effects. In Section 4 we apply our method to the Double 
Pulsar to obtain limits for preferred-frame effects from this 
system. Finally, wc present results of simulations that show 
how the precision of this test will improve within the next 
couple of years, and we investigate the possibility of combin- 
ing observations of several binary pulsar systems, to obtain 
a full sky coverage with high sensitivity. 



2 PREFERRED-FRAME EFFECTS IN THE 
MOTION OF COMPACT BINARIES 

In this first part wc study the motion of compact binary 
systems in theories of gravity that have a preferred frame of 
reference. We derive the equations of motion from a general 
formalism introduced by Will (1993) which includes strong 
field contributions using generic strong-field parameters. As 
we do not impose boost-invariance, the resulting equations 
of motion contain terms that depend on the motion of the 
binary system with respect to a preferred frame. Hence, our 
computations generalize the earlier work by Nordtvedt and 
Will (1972) and Damour and Esposito-Farosc (1992b) to de- 
scribe binary systems and their orbital parameters with ar- 
bitrary eccentricities. Using the method of perturbations of 
osculating orbital elements it is shown that the rate of rela- 
tivistic precession of periastron, the orbital eccentricity, and 
the orbital period should change in a characteristic way over 
time if a binary system moves relative to a preferred frame 
of reference. 



2.1 Binary dynamics in the generalised EIH 
formalism 

When dropping constant terms and rescaling the masses, 
the semi-conservative, generalised Einstein-Infeld-Hoffmann 
(EIH) Lagrangian of Will (1993) reads for a system consist- 
ing of two compact objects 



where 

r(0) 1 2,1 2 , 

' = -^rripVp + ^rUcVc + 
and 



QGrnpiric 



(1) 



(2) 



L^^^ = \ApmpVp + \AcmcV, 



2r 
Gm„mc 



[3B(u? + «p) - 7C(vp • Vc) - f (vp • ii)(vc • n)] 



2r2 



[rUpVc + nicVp] 



(3) 



rup and rric are the inertial masses of pulsar and companion, 
respectively, r = |r| = |xp — Xc| is the coordinate distance 
between pulsar and companion, n = r/r, Vp and Vc are the 
coordinate velocities of pulsar and companion, respectively, 
and G is the Newtonian constant of gravity. 

The coefficients Ap, Ac, B, C, Vp, Vc, £, and Q are 
functions of the parameters of the chosen theory of grav- 
ity and of the structure of each body. They account for 
contributions of the highly relativistic interior of the neu- 
tron stars to the binary dynamics. In particular, it is as- 
sumed that these coefficients are independent of the inter- 
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body distance in the binary system. In this paper we wiU de- 
note them as 'strong gravity coefficients'. The strong grav- 
ity coefficients B, C, £, and Q are symmetric under inter- 
change of p <-> c. For instance, in the fully conservative 



version of Rosen's bimctric theory Q — 1 



where 



Sp and Sc are the 'ffist sensitivities' (= —dlnm/dlnG ~ 
[gravitational binding energy] /[mass]) of pulsar and com- 
panion, respectively. In GR, these parameters are unity and, 
therefore, the compactness of a body does not have an im- 
pact on its orbital motion, a property of GR known as 
the 'effacement' of internal structure (see the discussion in 
Damour 1987). 

We note that combining the PFE- Terms of the La- 
grangian of Damour and Esposito-Farcsc (1992b) and 
the generalised conservative Lagrangian of Damour and 
Esposito-Farese (1992a) yields a Lagrangian that is equiva^ 
lent to the generalised Lagrangian of Will (1993) if one sets 
the ^-terms in Will's Lagrangian equal to in unity. As it 
turns out that these ^-terms play an important role in the 
violation of Lorentz invariance, we use Will's Lagrangian for 
our calculations. 

In contrast to Will (1993), we do not impose post- 
Galilean invariance when calculating the equations of mo- 
tion, in order to include gravitational theories that have a 
preferred frame of reference and, therefore, are not post- 
Galilean invariant. Hence, the equations of motion for pulsar 
and companion in the preferred frame are derived from the 
Euler-Lagrange equations 

d dL dL „ , d dL dL „ 

(4) 



dt dv„ dx„ 



• and — - — 
at avc 



= 0. 



without any further restrictions on the 'strong gravity co- 



efficients' of L. The relative acceleration a 



a.c for a 



binary pulsar system at rest with respect to the preferred 
frame is then found to be 



a = a(°)+a« 



where 



(0) fn n 
a = TT- 



I 2 



(5) 



(6) 



m 
r 



-(i-©:)^p-(i-D;)e 

- i(3r-l) + i(6 + at + 02)5^^0 

+ |(l + «2)6^c(v.n)^} 
+ ^ [l + 3B* - (2 - at + aV)ipi, 

-(i-A)e^^(i-^p)dl (v-ft) 



(7) 



where m* = QG{mp+mc), S,p = mp/m, = nic/m = 1— ^p, 

B* = B/g, c* = c/g, v; = Vp/g^ v* = v^/g^ s* = s/g, 

and further 



al=£* + 7C* -QB* -2, 
a*2= £* - \ . 



(8) 
(9) 



Note, in GR a\ = ot2 = 0. Further, a{ and a2 are propor- 
tional to the PFE parameters qi and 0.2 of Damour and 
Esposito-Farese (1992b), respectively: a* = g&i. 

The 'Newtonian' part of the above equations of motion, 
a^°^ , has the well-known Keplerian solution 



P 



r = r(ex cos (^-t-ey sine/)) , r = 

1 -h e cos { 

where 

p = a{l — e^) and (Pj, /27r) ^ = a^/m* . 



(10) 

(11) 

(12) 



Pb is the orbital period of the binary system and e the ec- 
centricity of its orbit. 

The post-Newtonian terms, a'^-'/c^, produce a secular 
advance of periastron given by 



(13) 



where 

p = r + i(i-D;)Cc + i(i-o:)ep 

+ |(2at - q^)CpCc - i(l - Ap)C - i(l - AcKp . (14) 



2.2 Preferred-frame effects in the relativistic 
binary motion 

In the previous Section we have given results for a binary- 
pulsar system that is at rest with respect to the preferred 
frame of reference. In this Section we investigate the dynam- 
ics of a binary system that moves relative to the preferred 
frame with velocity w. As a matter of convenience, we will 
use a frame that is comoving with the binary system. The 
equations of motion in the comoving frame can be derived 
by imposing a post-Galilean transformation (Chandrasekhar 
and Contopoulos 1967) on the equations of motion that re- 
sult from the Euler-Lagrange equations (4) of the previous 
Section. For the relative acceleration, as expressed in the 
comoving frame, one finds 



JO) 



+ a«/c'+a^™Vc 



(™) / 2 



(15) 



where 



(w) 

a^ ' = 



2r2 



[«i(Cp - Cc)(w ■ v) - 3a; (w • fi)^ 

+ {1 - Ac)£.p'w^ + (1 - Ap)^cw'^ 

-(1 - Ac)ep{^ • v) + (1 - A)d(w • v)] 

+ [Qt(Cp - ?c)(v • n) - 2a; (w • n) 

-2(1 - ^c)Cp(w • n) - 2(1 - Ap)ic{w ■ n) 
+2(1 - A)ep(v • A) - 2(1 - A)d(v • &)] 

[(1 - • fi) - (1 - Ap)e,{w . n)] , (16) 

The variations of the orbital parameters of a binary pul- 
sar caused by the preferred- frame acceleration a'-™'/c^ can 
now be calculated by the standard technique of perturba- 
tions of osculating orbital elements. Using the notation of 
Damour and Taylor (1992), one finds, when averaging over 
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one full orbit, for the change in the longitude of periastron 



Aw'"'' = 



smx 



sin ?/)cos2x 



■Ancosi, (17) 



for the change in the longitude of the ascending node 



AO*™' 



c^p sin I 

2 



Qie — cos -i/) cos a; 

"0/ 



-Q2 ( — ) sin2ip 

\VoJ 



cos X sm a; 



+ sin X cos a; 



for the change in the projected semi-major axis 
Ax*""' Trm 



(18) 



X c^p tan I 

+0.2 (^)%in2V' 



Qie ( — ) COST'S 
cos X cos w 



sm X sin w 



(19) 



for the change in the eccentricity 



AeW = ^(l + eVe 



Qi ( — ) sin-!/)cosx 



+Q2 



a - e2 Vvo/ 



sin Vsin2x 



(20) 



for the change in the mean anomaly (dropping terms inde- 
pendent of x) 



AM' 



Aw*'"' -h An'"'' cos? 



12QieFei/l-e2 (—\ sin V sin x 
+Q2 (2FeV'l-e2- 1) sin^^cos2x 

where 

Qi = al (^p - ^c) - 2(1 - Aj,)e, + 2(1 - Ac)fp , 

Q2 = 02 + (1 - Ac)^p + {1- Ap)ic , 

Q'l = (1 - Ap)e, - (1 - A)ep , 

Q'2 = a;- 2(1 - - 2(1 - ^p)ec , 

and 

1 



Fe = 



(21) 

(22) 
(23) 
(24) 
(25) 

(26) 



■0 is the angle between k, the direction of the orbital angular 
momentum, and w. x is the angle between the periastron of 
the pulsar and the projection of w into the orbital plane. 

In fully conservative theories of gravity one finds Ap = 
Ac = 1 and Ui = 0.2 = 0. Consequently fully conservative 
theories of gravity do not predict any preferred-frame effects 
in the motion of compact binaries, i.e. Qi = Q2 = Q'l = 
Q2 = 0. In GR, in addition to the absence of preferred- 
frame effects, one has V = 1. 

2.3 The PPN limit 

In the true post-Newtonian limit for masses with negligible 
self gravity one finds the strong gravity coefficients as func- 
tions of the PPN parameters (3, 7, ai, and 02 (Will 1993): 



G = l, 

r = 1(27 + 1), 

C* = i(47-|-3-|-ai-a2), 

v; = vi = 2i3-i, 

£* = 1 + a2 , 
and 

al = ai , 

0:2 = Q!2 ■ 

Hence 

P = i(2 + 27 - /?) + i(2ai - a2)^p^c , 
Qi = ai(Cp -Cc) , 
Q2 = a2 , 
Q'l = . 
Q2 = 0:2 , 



(27) 
(28) 
(29) 
(30) 
(31) 
(32) 

(33) 
(34) 

(35) 
(36) 
(37) 
(38) 
(39) 



Nordtvedt and Will (1972) derive the contribution of a 
preferred frame to the advance of periastron and the change 
in the eccentricity for planetary orbits (equations 61, 62). In 
the PPN- limit our equations (17) and (20) agree with their 
results if one assumes nip <^ rUc and e 1. As we restrict 
our discussion to semi-conservative theories of gravity, we 
do not get the self-acceleration terms related to 03. 



3 TESTING PREFERRED-FRAME EFFECTS 
IN BINARY PULSARS 

Wo have seen in the previous Section that if a binary pulsar 
moves with respect to a preferred frame of reference, changes 
of its orbital parameters occur, which should then become 
apparent in its timing data. Therefore, assuming a direction 

and velocity w for the preferred framed we will demonstrate 
that timing data can be used to constrain or even measure 
the strong field parameters Qi and Q2, which determine the 
strength of the preferred-frame effects. 



3.1 Measuring PFE efTects 

Equations (17) to (21) can be converted into (averaged) 
first order time derivatives in the orbital elements via (w) = 
Auj/Pb, etc. At a first glance, comparing these expressions 
with the observed values for cu, x, and e ^ obtained in tim- 
ing observations, seems to provide a method to measure or 
at least constrain the strong- field parameters Qi and Q2. 
In particular, using a measurement for u! seems to be very 

^ We do not make any assumption about the nature of this pre- 
ferred frame in our derivation. Only later, we investigate obser- 
vational data assuming some natural choices such as the motion 
of the binary pulsar with respect to the cosmic microwave back- 
ground, representing a distinguishable frame of reference in the 
Universe (Nordtvedt and Will 1972). 
i Note that CI is not an observable quantity. 
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attractive, as this timing parameter is usually the easiest 
to determine among the so called "Post-Keplerian" param- 
eters. However, this seemingly straightforward approach of 
using only the first order time derivatives has a number of 
problems: (a) the gravitational mass parameter m* can usu- 
ally not be determined in a theory-independent way, (b) the 
value of the strong-field parameter V is not known, hence 
prevents the use of (cLi) , and (c) for the most promising test 
systems the orbit is usually relativistic enough to cause a 
non-linear evolution of the orbital parameters. We discuss 
each these potential problems in turn. 

Without making detailed assumptions on the strong- 
field properties of the underlying theory of gravity, most bi- 
nary pulsars will not allow the determination of m* and the 
orbital inclination i as required in equations (17) to (21). 
Indeed, most Post-Keplerian parameters, like the Einstein 
delay, 7, or the change in orbital period due to gravitational 
wave damping, Pi,, cannot be used as they are expected 
to have strong-field contributions, which are not quantified 
within a theory-independent framework as used in this pa- 
per. In contrast, under some very natural assumptions con- 
cerning the leading terms in the space-time metric in the 
inter-body region and far from a binary pulsar system, the 
observation of a Shapiro delay can be used to determine the 
sine of the orbital inclination. We will discuss this point in 
more detail later, when we also demonstrate that under very 
special circumstances, such information can be used to ac- 
tually determine both mass parameters m* = QGrUp and 
ml = QGrric. 

Even if we could determine m* from post-Keplerian 
parameters not involving o), we still could not extract the 
preferred-frame contribution in {w) as we do not have 
a (theory-independent) value for V of equation (13). In 
other words, usually one cannot separate the preferred-frame 
terms from the other relativistic terms. 

Moreover, one can precisely measure post-Keplerian pa- 
rameters only in relativistic binary pulsars. In such a case, 
the preferred- frame efi'ccts as described by equations (17) 
to (21) cannot be described by a simple lincar-in-time ex- 
pressions, like e = {e){t — to). Fitting for first order time 
derivatives, therefore, cannot be used to consistently test 
for preferred-frame effects in these systems. This is true for 
any binary pulsar system where due to a relativistic ad- 
vance of periastron the longitude of periastron has changed 
significantly since the time of its discovery. For instance, in 
PSR B1913-fl6 u! has advanced by about 135 degrees since 
its discovery in 1974, and therefore first order time deriva- 
tives cannot be used to constrain preferred frame effects, as 
done in Bailey and Kostelecky (2006). The result is a signif- 
icant change in x since x = const. — Lij{t). Consequently, the 
time derivatives of the orbital parameters are trigonomet- 
ric functions of a;(f) and, accordingly, the preferred-frame 
effects impose periodic changes on the orbital parameters 
themselves, with frequencies u> and 2uj. 

In order to test for preferred-frame effects in such sys- 
tems, one needs to include this periodic 'signature' of a pre- 
ferred frame into the timing model, which is used for fitting 
the timing data. In the following we will describe a timing 
model and apply it to the Double Pulsar. 



3.2 A timing model for the 'signature' of 
preferred-frame effects 

3.2.1 Case 7. i ~ 90 degrees 

As outlined in the previous Section, the presence of 
preferred-frame effects in the motion of a binary system re- 
sults in periodic changes of the orbital elements. In this Sec- 
tion we will extend the existing standard Damour & Dcru- 
ellc (DD) timing model (Damour & Deruelle 1985, 1986) to 
include these periodic changes and, therefore, allows to fit 
for the amplitudes of these preferred-frame effects. We will 
demonstrate our method by applying it to the Double Pul- 
sar system, and so we restrict the following discussion to 
binary systems with an orbital inclination, i, close to 90 de- 
grees. A generalisation for systems with orbital inclinations 
significantly less than 90 degrees is straightforward. 

If « ~ 90° changes in the projected semi-major axis 
(equation 19) arc small and can be neglected in the timing 
model. In the following xo denotes the angle between the as- 
cending node and the projection of w into the orbital plane, 
wx, and consequently x = Xo — <^ (see Fig. 1). 

Integrating equation (17) for i = 90° one finds 



aJ^^ (T) = nl'^' cosiuL - xo) - V^^' sin 2{u;l - xo) , (40) 



where 

wi, =a;o+w^'^(r-ro) , 
and 



,M = ^(i-eF.^) (™)sin^, 
12'P Ve J \voJ 

fwY . 2 , 
— sin tp . 

\VoJ 



^ 12P " \vo 



(41) 



(42) 



(43) 



When integrating, we have made the assumption that it is 
sufficient to keep terms linear in r|^"\^ Hence uj^^^ is given 
by equation (13). 

Integrating equation (21) one finds 

Ae*"'' (T) = ni'^ sin{^L - Xo) + vi"^ cos 2(0;^ - xo) . (44) 
where 

Si 



= + - sinV>. 

(e) _ Q2 e(l-|-e^)j; fw^ 2 
V2 — 



(wy . 2 , 

— sm w . 
\vo/ 



(45) 
(46) 



1273 ^/T^ 
Integrating equation (20) one finds 

AM*'"' (T) = cos(a;i - xc) + vi^^ sin 2{ujl - Xo) ,(47) 

where 



(M) (lj) n z , 



9± 
V 



Jl-e'i — VsinV, (48) 
\voJ 



(M) 
»72 



Q2 



(^eV^-|)( 



sin'^V ■ 



(49) 



§ It can be verified directly by fitting the timing data if this 
assumption is applicable. 
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V 

Figure 1. Definition of angles in tiie PFE timing model. 

We emphasize that the ratios ??i"V'?i'^^ ^'^^ vi"'^/V2'^ are 
only functions of the Keplerian eccentricity, e, of the binary 
system. 

For a consistent, theory- independent analysis of binary 
pulsar timing data, Damour and Deruelle developed a phe- 
nomenological timing model, using a parametrized post- 
Keplerian approach (Damour & Deruelle 1985, 1986). The 
DD timing model applies to the class of Lorcntz-invariant 
theories of gravity, such as the well-motivated tensor-scalar 
theories. However, it does not describe the preferred-frame 
eflfects that we consider here. On the other hand, it is easy 
to incorporate these variations into the DD timing model by 
simply adding solutions (40), (44), (47) accordingly to the 
the DD model: 

uj{T) ^ Lu{T) + Alu{T)'-'"\ (50) 
e(T)^e(T) + Ae(r)("', (51) 

^{T-To)^^{T-To) + AM{tY'^K (52) 

Therefore our new timing formula, for binary pulsars with 
i ~ 90°, contains, in addition to the Keplerian and post- 
Keplerian parameters of the DD timing model, five new 
preferred-frame timing parameters: the four PFE ampli- 
tudes rii"\ vi^K Vi^\ Vt^K and the angle xo- The ampli- 
tudes ri[''^ and ri2^^ are proportional to ri["-' and rjj'^-' respec- 
tively, where the factors are only functions of the Keplerian 
eccentricity, e. 

If the longitude of periastron has advanced for a con- 
siderable amount (^ 90°) since the discovery of a binary 
pulsar, the timing data will allow to fit for all of these new 
timing parameters. If the binary motion shows the effect 
of a preferred-frame in the Universe, all five independent 
'preferred-frame parameters', ri["\ Vi"^ 'n2^\ and xo 

can be determined. One can even decide to fit additionally 
for ril"^ and r^j*' i-e. perform a fit simultaneously for these pa- 
rameters and for ri['^^ and 772"'' • The existence of a preferred 
frame could then be confirmed by verifying the expected 
relationship between the PFE amplitudes in uj and e. The 



ratios of these amplitudes are unique in the sense that they 
are a function of e only and therefore independent of any 
strong-field parameters, and of the magnitude and direction 
of w. For the parameters of the Double Pulsax we would 
expect to measure 

vi'^^/vi'^ = 11.262 , (53) 

r?('^V^2'' = 5.642 (54) 

if a preferred frame exists. This first tost would incontro- 
vertibly reveal the existence of a preferred frame. We will 
show later, how the direction of the preferred can also be 
constrained or even determined. 

In the absence of any preferred- frame effect, we expect 
none of the PFE timing paramotors to have significant val- 
ues. In particular, the magnitude of PFE amplitudes should 
be consistent with zero, and the angle xo is undefined. In this 
case, one can obtain upper and lower limits for ti["^\ V-z^K 
r/j*^' , and r/j^^' by iiolding xo fixed to a particular value 
while fitting for the PFE amplitudes. Stepping through the 
possible values xo £ [0, 27r], various directions on the sky 
can be probed. In the case of a relativistic binary pulsar, 
a clear separation of the preferred- frame contribution (r/J"' , 
jjj'^') from the other relativistic terms (a''^') in the preces- 
sion of periastron is possible and will improve with time due 
to increasing coverage of w-space. 

3.2.2 Case II. i ^ 90 degrees 

As stated earlier, the above assumes that i ~ 90°. If one al- 
lows for any orbital inclinations, one finds that equation (40) 

is replaced by 

Ao; = cos{LOL-Xo + 5'i)-r]i,'^'' sin 2(u;i - xo + <52) .(55) 

Like ?7i"' and ri2'^\ and 773"' are proportional to Qi/V 
and 0.2 /V respectively, while 5^ and 82 depend only on Ke- 
plerian parameters and the orientation of the pulsar with 
respect to the preferred frame. In this case also x{T) has to 
be considered, for which one finds an expression similar to 
equation (55). Following the procedure above, a derivation 
of these expressions is straight forward. 

3.3 Translating timing model parameters to 
strong-field PFE parameters 

In this Section we give a description on how to convert the 
timing results, obtained after applying the new PFE timing 
model., into the strong-field PFE parameters introduced in 
Section 2. 

As discussed in the previous Section, for a binary pul- 
sar with i ~ 90° values for the timing parameters ri^\ 
ri^\ ri[^''\ and 773*^^ can be converted into values for 
Qiw sin wsinip/V, Q^w^ siv?' ip /V . 

Without any further assumptions, however, we cannot de- 
termine or even restrict the angle ip. Here one has to keep in 
mind that, in general, the longitude of the ascending node of 
the pulsar orbit, O, cannot be determined from pulsar tim- 
ing observations and has to be treated as a free parameter. 
We discuss possible exceptions further below. 

Relating the orientation of the orbit to the movement 
relative to a preferred frame in direction w, requires the 
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knowledge of the orbital inclination angle. In cases where the 
binary orbit is seen nearly edge-on, one can often determine 
this angle accurately, as one expects to observe an additional 
delay in the arrival time of the pulsar signal during superior 
conjunction of the pulsar, caused by the gravitational field 
of the companion. This is the case in the Double Pulsar 
where this "Shapiro effect" cam be used to determine the 
orbital inclination i (modulo the ambiguity i ^ tt — i). In 
order to do this, one only have to make some general, very 
natural assumptions. As in Will (1993), we assume that to 
first order the metric valid in the inter-body region and far 
from the system can be written as 



goo = -1 + 



goj 



9ij = 1 + 



2GK*mr, 



+ ■ 



(56) 
(57) 
(58) 



where k* , 7* , k* , and 7* are functions of the parameters of 
the theory and of the structure of pulsar and companion. In 
GR = 7p = «* = 7* = 1. The time taken by the pulsar 
signal to travel from the pulsar to the solax system in such 
a metric can be accounted for by adding 



■In 



e cos U 



sin i j^sin c<;(cos U — e) 
+ (l-e^)^''^cosa)sin?7 | . (59) 



to the timing model, which can be used to fit for sin i, the 
"shape" of the Shapiro delay, without any knowledge of 
the strong-field parameters (Blandford and Teukolsky 1976, 
Damour and Deruelle 1986b, Will 1993). U is the eccentric 
anomaly as defined by the DD model. The measurement of 
sini then allows to exclude certain directions of w in the 
sky, given certain values of xo- In fact, if one makes certain 
assumptions of the magnitude and direction of w, i.e. the 
direction and motion of the preferred frame, a measurement 
of sini (and the computation of xo) will allow the deter- 
mination of sintp, and, consequently, the determination of 
Qi/P and Qi/V for the assumed preferred frame from the 
fitted PFE amplitudes. 

On the other hand, as outlined earlier, if no preferred- 
frame effects are observable in the orbital motion of a bi- 
nary pulsar, fitting for the PFE amplitudes (while hold- 
ing Xo fixed) can be used to determine limits on Qiw/V, 
Q.ivP' jV , Q!\w/V, and Q'^w^ /V for (nearly) any direction 
in the sky. However, during this transformation from PFE 
timing model parameters to these physical quantities, one 
needs to determine xo and t/j for a given direction in the sky. 
That requires the knowledge of the, in general, unknown an- 
gle n, while one also has to account for the ambiguity in the 
sense of the inclination i. Two practical approaches exist to 
overcome this problem. These angles can either be chosen 
such that the limits are most conservative, or one can per- 
form Monte-Carlo simulations varying over Q, and the two 
possible values for i. Assuming a uniform distribution for Q, 
and equal probability for i and n — i, it is possible to deter- 
mine limits for any chosen confidence level. It is clear that 
the most conservative method cannot provide any restric- 
tions for directions along a cone that has an opening angle 



i around the line of sight, as for any direction on this cone, 
Q can be chosen such that 'ip = 0. 

The situation is obviously improved if can be deter- 
mined. We note that the proper motion of a binary pulsar 
can produce secular changes in the orbital elements that de- 
pend on the longitude of the ascending node, f2, and the 
orbital inclination, i, as it changes the apparent geometri- 
cal orientation of the orbit (Arzoumanian et al. 1996 and 
Kopeikin 1996). In principle these changes can be used to 
determine and i from timing observations. However, in 
practice it will be difficult to separate the proper motion 
effects from the relativistic changes of the binary orbit, in 
particular if we make only very generic assumptions about 
the underlying theory of gravity, as we do in this paper. 

In principle, measurements of the timescale of the in- 
terstellar scintillation (ISS) over an orbit can also be used to 
estimate the orbital inclination i and the transverse velocity 
of the centre of mass of the system (Lyne & Smith 1982, 
Ransom et al. 2004, Coles et al. 2005). A comparison of the 
transverse velocity derived from ISS and the transverse ve- 
locity obtained from timing observation could then be used 
to determine Q. However, the scintillation-based velocity de- 
pends on a number of assumptions about the properties of 
the effective scattering screen and, therefore, are more sus- 
ceptible to systematic errors than timing measurements. 



4 THE DOUBLE PULSAR 

Following the discovery of the Double Pulsar system in April 
2003 (Burgay et al. 2003; Lyne et al. 2004), timing obser- 
vations have enabled the most stringent tests of CR in the 
strong-field regime to date (Kramer et al. 2006). Indeed, the 
system is a unique laboratory for gravitational physics for 
a number of reasons. Firstly, both members of the binary 
system are visible as active radio pulsars, providing access 
to a simple measurement of the ratio of their masses and 
hence providing theory-independent (at least to IPN order) 
constraints for tests of theories of gravity. 

Secondly, the fortunate orientation of the binary orbit in 
space allows us to observe the Double Pulsar under a nearly 
perfect edge-on geometry. In addition to precise measure- 
ments of a Shapiro delay, this enables further independent 
estimates of the system's inclination angle. ^ 

The combination of the determination of sini via the 
Shapiro delay with the measurement of the mass ratio al- 
lows a theory-independent determination of both mass pa- 
rameters, m* and m*, as we will discuss in detail later. In 
this case V can now be measured independently, thus di- 
rectly giving limits for Qi and Q2 from fitting for J7j"^^ and 
ri2^^ . Further, from this one can determine al and 02 for 
gravity theories that fulfil Ap = Ac = 1. In this sense the 
Double Pulsar is indeed a unique laboratory for strong-field 
preferred-frame effects. 



^ Observations and modelling of the eclipse of the radio emis- 
sion of the millisecond pulsar PSR J0737— 3039A during superior 
conjunction, lasting for about 25-30 seconds, provides another es- 
timate (Breton et al. 2006), while two methods using observed 
scintillation properties provide two others (Ransom et al. 2004; 
Coles et al. 2005). 
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Relative angle Xq C™*^) 

Figure 2. Limits on the existence of preferred-frame effects 
present in the orbital motion of the Double Pulsar as measured 
by the parameters Tj^"' (dark grey) and r/j"^' (light grey) as in- 
troduced in the PFE timing model. See text for details. 



The third reason why the Double Pulsar is a unique 
gravity lab is the compactness of the orbit. With larger mean 
orbital velocities than in any other binary pulsar, the system 
is the most relativistic system known. This is reflected by the 
measurement of the largest advance of periastron observed 
in any binary system. As a result, the Double Pulsar orbit 
precesses by 16.9 degrees per year with respect to any possi- 
ble existing preferred frame, changing the angle x in Fig. 1 
by a corresponding amount. As such the Double Pulsar is an 
ideal candidate to study the existence of preferred-frame ef- 
fects in strong-field gravity as the magnitude of co does not 
only allow its very precise measurement (currently at the 
4 X 10~* level) but the orbit also 'covers' a lot of periastron 
angle space in a rather short time span. According to the 
theoretical framework laid out in the previous sections, such 
a system is ideal for testing for PFE variations in the orbital 
parameters, the 'signature' of a preferred frame. 



4.1 Application of the PFE timing model 

In Section 3, we have developed the PFE timing model that 
introduces, based on the DD timing model, the new PFE 
model parameters ri['^\ Vi^'K t}^''^ Xo- We now want 
to apply this model to the timing data of the Double Pulsar. 
Wo make use of nearly three years of timing data obtained 
by Kramer et al. (2006) during which the pulsar periastron 
has advanced by about 50 degrees. 

At present the timing data of the Double Pulsar system 
does not show a PFE 'signature'. Hence, we can use the data 
only to derive limits for the PFE amplitudes by fitting for 
the PFE amplitudes while holding xo fixed at a given value. 
We then vary xo in sufliciently small steps between and 
360° to derive limits for any value of xo- Fig. 2 illustrates 
the 95% confidence for rf^' and ri"^\ consistent with the 
non-existence of preferred frames. 

In principle, we would also aim to fit for 77j*^^ and 7]^^ 
which would yield limits for Q'l and Q21 giving separate 
limits for a.\, a.2, Ap, Ac when combined with the results for 



77}'^ and ■ However, the Double Pulsar timing data used 
in this paper do not yet allow to fit for rif[^^ and r)^^ . The 
periastron has not advanced far enough to separate these 
PFE amplitudes from other orbital parameters, in particular 
the decrease in the orbital period due to the emission of 
gravitational waves. As we will show later, with longer time- 
spans and more coverage of the periastron angle, such a 
fit will be possible, while simultaneously the amplitude of 
all limits will be greatly reduced. It is obvious that some 
orientations of the orbit with respect to w lead to better 
constraints than others. 

4.2 Determination of quantitative limits 

In this Section we translate the limits on the 'signature' of 
a preferred frame into quantitative limits within the gener- 
alised EIH formalism. Our goal is, first, to determine limits 
for QiU! and Q.2W^ for any direction in the sky, and secondly, 
to determine limits for Q\ and Q2 for a preferred frame that 
is at rest with respect to the CMB and one that is at rest 
with respect to our galaxy. 

The first step is to determine the masses from the timing 
data described by applying the PFE timing model. Equa- 
tion (59) can be used to convert the measurement of the 
Shapiro "shape" timing parameter into a measurement of 
sini. Now, combining equation (3.15) of Damour and Tay- 
lor (1992) II with the measurement of the mass ratio obtained 
from R in the Double Pulsar, we obtain the following values 

nil = QGrUp = (1.339 ± 0.002) x GMq , (60) 

m* = GGrric = (1.250 ± 0.002) x GMq . (61) 

It is at this point where we exploit the uniqueness of 
the Double Pulsax which provides us with measurements 
of m* and rn*. Combining these with the observed rate of 
periastron advance allows us to determine the strong-field 
parameter V (sec equation 13). The result is 

P = 1.000 ±0.001 . (62) 

In order to establish the absolute orientation of the bi- 
nary system in space, one needs the orientation of the as- 
cending node n, which is still unknown for the Double Pul- 
sar (but see discussion in Section 3). As described earlier, we 
perform Monte-Carlo simulations and vary this parameter 
between and 2n for any given direction of w. For a given 
w and a given fi, we compute the angle xo and read off the 
corresponding limits for r][^^ and TJa"^^ from the calculated 
and tabulated values shown in Fig. 2. To account also for the 
ambiguity in the sense of the inclination i, the whole proce- 
dure is performed for 10,000 datasets for each direction in 
the sky and inclination angles of i and tt — i. As discussed 
earlier, the obtained values for r]['^^ and rji^^ can be trans- 
lated into limits for Qiw and Q.2w'^ using equations (42) and 
(43). Figs. 3 give the results of our numerical simulations for 
preferred frames associated with different directions of the 
sky as seen from the Double Pulsar. The limits shown as 
a colour map are derived for an assumed binary pulsar ve- 
locity relative to the preferred frame of 100 km/s. Results 

II Note that the effective gravitational constant Q in Damour and 
Taylor (1992) corresponds to QG in this paper. 
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Figure 3. Limits computed for 95% confidence levels on lQi| and IQ2I for the sky as seen by the Double Pulsar. The limits are calculated 
for |w| = 100 km/s, and have to be scaled accordingly for other values of |w|. Directions discussed in the text are marked by the blue 
area, i.e. direction relative to a Galactic frame (top left) and direction relative to the Cosmic Microwave Background (CMB) (right). The 
extension of these areas reflects the uncetainty in the radial velocity of the binary pulsar system with respect to the solar system. When 
extracting the limits for the Galactic and the CMB frame from these figures, one needs to keep in mind that the |w| of these frames 
changes across these areas. The location of the Sun as seen from the Double Pulsar is marked by 0. 
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Figure 4. Future improvement of PFE tests in the Double Pulsar. (Left) Amplitude of the PFE parameters (dark grey) and jjj'^' 
(light grey) as expected for 10 years of Double Pulsar timing data. (Right) Evolution of limits on »7j"' and r/j^^ as a function of time by 
recording their worst values over all xo at a given epoch. The year of observation is indicated at the top of the figure. 



for other relative speeds can be scaled according to equa- 
tions (42) and (43). 

We also indicate the directions for two selected frames 
that may be considered as being related to a preferred frame; 
we indicate the direction of motion relative to the CMB and 
the direction of the motion relative to the Galactic reference 
frame. The motion of the solar system with respect to the 
CMB we take from Hinshaw et al. (2006) and Lineweaver et 
al. (2006), and with respect to the Galactic reference frame 
from Mignard (2000). To calculate the motion of the Double 
Pulsar with respect to these frames of reference we need to 
know the motion of the Double Pulsar with respect to the 
solar system. While a proper motion of the pulsar has been 
measured, one of the few unknown system parameters is the 



undetermined radial velocity of the system. In accordance 
with recent evolution studies of the Double Pulsar system 
(Stairs et al. 2006), we perform our analysis for a range of ra- 
dial velocities between —100 and +100 km/s, and therefore, 
the direction relative to the CMB and the Galactic frame 
convert into small areas on the sky, from which we pick the 
worst 95% C.L. limits on Qi and Q2: 



CMB: -0.02 < Qi < 0.01 
GAL: -0.01 < Qi < 0.01 



-0.3 < Q2 < 0.2 
-0.3 < Q2 < 0.3 



(63) 
(64) 



For gravity theories with Ap = Ac = 1 we obtain directly 
limits for al and 0:2 (see equations 22 and 23): 



CMB: 



-0.5 <al < 0.3 -0.3 < al < 0.2 



(65) 
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GAL: -0.3 <al< 0.3 -0.3 < < 0.3 (66) 

The limits given here arc not very tight yet. In fact, 
they are several orders of magnitude worse than the lim- 
its for the weak-field PPN parameters ai and 02- On one 
hand, as we will discuss below, these limits will improve 
considerably during the next couple of years. On the other 
hand, these limits hold for preferred-frame effects related to 
strong gravitational fields. If one expands, for instance, a* 
as a function of the sensitivity of the gravitation body 

al = ai + Qii {sp + Sc) + ... (67) 

one sees that tests in the solar system made to restrict ai, 
can only test terms which are not related to the internal 
gravitational fields of a body. Testing those strong-field co- 
efficients, like Qii, has to be done via binary pulsars. 

Our strong-field limit for al has to be compared with 
the limits obtained by Damour and Esposito-Farese (1992b) 
and Wex (2000), which are better by about three orders of 
magnitude. However, in a few years from now wc expect the 
limits obtained from the Double Pulsar to be comparable 
with the limits obtained from small-eccentricity binary pul- 
sars, and as outlined in Section 3 the method presented here 
does not need to rely on probabilistic considerations. 

Emphasizing the arguments and results presented by 
Will (1993) and Damour and Esposito-Farese (1996a), we 
point out that Double Neutron star Systems (DNSs) are par- 
ticularly sensitive to theories of gravity where higher order 
terms in Eqn. (67) become important, so that for instance 
al oc SpSc (similarly, 02 oc SpSc), since for a white dwarf 
SWD ^ 10~^ < sns. 

4.3 EYiture measurements 

The continuing precession of periastron in the Double Pul- 
sar will lead to a constant improvement in the limits on pre- 
ferred frames using our proposed method. In order to predict 
the future limits we have simulated timing data for the Dou- 
ble Pulsar, assuming similar timing precision and frequen- 
cies of observations as described in Kramer et al. (2006). 
Simulating data until 2013, hence for 10 years after the dis- 
covery of the system and half of the periastron precession 
period, we have computed the limits on the observed param- 
eters ri["^ and r^j"^ as a function of time. In these simulations 
we make the unlikely and hence conservative assumption 
that the timing precision will fail to improve over the next 
years, so that the quoted limits should be considered as con- 
servative also. The result is shown in Fig. 4. Both parameters 
scale essentially similarly, and in 2013 the limits will have 
improved by more than two orders of magnitude, resulting 
in limits of |qS;| < 2 x 10"^ and loa] < 2 x 10"^. The vast 
improvement is also demonstrated by comparing Figure 2 
with the left part of Figure 4 where we plot the limits on 
77'"' and 772"^ as a function of xo as expected for the year 
2013. These predicted values are still somewhat worse than 
those derived by Damour and Esposito-Farese (1992b) and 
Wex (2000) , but we estimate that with the timing precision 
being certain to improve over the next few years (e.g. by 
using the telescopes to be constructed as pathfinders to the 
Square-Kilometre Array (SKA)), the limits will be improved 
upon further. 

In addition, by 2013 the periastron will have advanced 



by nearly n since the pulsar discovery, which should also 
allow a separate measurement of rj['^''^ and r;2*^^ with com- 
parable precision. This will provide us with measurements 
for Qi, Q2, Q'l, and Q'2, and consequently with values for 
the strong- field parameters at, 02, Ap, and Ac. 



5 AN PFE ANTENNA ARRAY 

In the future, further binary systems with a variety of ori- 
entations in space will be discovered. In particular, with the 
Square-Kilometre Array (SKA) we expect to find about 100 
relativistic DNSs (Cordes et al. 2004; Kramer ct al. 2004). 
Some of these new systems will also be suitable sources for 
tests of local Lorentz-invariance of gravity, so that their ob- 
servations can be combined to derive further constraints, as 
also pointed out by Bailey and Kostelecky (2006). Indeed, 
the sensitivity of additional systems to difFcrcnt directions 
in the sky will difi'er from that of the Double Pulsar, so that 
we can construct a dense network of antennae for studies 
of preferred-frame effects. This results in a smooth and high 
sensitivity toward preferred-frame effects over the whole sky. 

In order to demonstrate this idea, we have sinmlated 
timing data for the Double Pulsar system if it were located 
at a different position in the sky, namely the Galactic Cen- 
tre. Such a system would be sensitive to directions which 
complement those of the real Double Pulsar. This can be 
seen in Figure 5 where wo show this situation simulated for 
the timing data as expected in 2013. Comparing the left 
panel of this figure with Figure 3 again demonstrates the 
vast improvement in the PFE limits over the next five to 
seven years. The impact of the artificial system's difi^crent 
position in the sky is clear from a comparison of the left and 
middle panel. A combination of the two sets of limits for 
any given direction in the sky, as shown in the right panel, 
demonstrates the concept of a PFE antenna array. 

Wc point out that for tests of the most general the- 
ories only data from systems with a similar combination 
of pulsar and companion mass can be combined. This be- 
comes clear when considering that the parameters Qi and 
Q2 contain an explicit and implicit dependence on the neu- 
tron star masses. The implicit dependence arises from the 
fact that the involved strong-field coefficients are functions 
of the compactness of the bodies and hence of the masses 
and the equation-of-state of the neutron stars. 

If preferred-frame effects are indeed present in the tim- 
ing data of a number of binary pulsars, this array of pulsars 
could then be used to significantly restrict or even determine 
the direction of the preferred frame (modulo tt related to a 
change of the sign in the measured PFE amplitudes), as a 
measurement of xo in one such system excludes the direc- 
tion to lie within a cone with opening angle n/2 — xo around 
the line-of-sight to the pulsar. 



6 SUMMARY 

We have developed a consistent methodology to measure 
preferred-frame effects (PFE) related to the strong internal 
gravitational fields in relativistic binary pulsars. We made 
only very general assumptions about the underlying theory 
of gravity by using the semi-conservative generalised EIH 
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Figure 5. Limits computed for 95% confidence Hmits on |Qi| and IQ2I for the sky as seen by the Double Pulsar in 2013 (left) and a 
simulated Double Pulsar-like system located in the Galactic Centre (middle). Combining the independent constraints provided by the 
two systems (right) makes this PFE antenna array sensitive to most directions in the sky. For simplicity we have assumed that the two 
systems have a negligible relative velocity when combining the two data sets. It is clear that the addition of further systems will lead to 
a smooth coverage of the whole sky. 



formalism of Will, which incorporates strong field effects in 
the post-Newtonian motion of binary pulsars. We were able 
to show that in relativistic binary systems with a high rate 
of advance of periastron, a preferred frame will cause distinc- 
tive periodic changes in the orbital elements of the binary 
system. 

The newly developed PFE timing model extends the 
DD timing model by including these periodic changes, and 
thus having the amplitudes of these changes as additional 
timing parameters. We described in detail how the measure- 
ment of such PFE amplitudes can be converted into mea- 
surements of the strong field parameters related to preferred- 
frame effects. 

For inclination angles i ~ 90°, the PFE amplitudes for 
changes in the longitude of periastron and the eccentricity 
are related by a factor which is independent of strong field 
parameters and parameters related to the preferred frame. 
This can be seen as a unique "fingerprint" of a preferred 
frame. 

We have also shown that in the absence of preferred- 
frame effects, our formalism can be used to determine limits 
for these amplitudes, by this restricting the strong field pa- 
rameters related to the violation of the Lorentz invariance 
for nearly any direction in the sky. We demonstrated that 
the Double Pulsar is the ideal test system for preferred-frame 
effects in strong gravitational fields, and have presented first 
preliminary results, which however are clearly less stringent 
then present limits from small-eccentricity binary pulsars. 
On the other hand, simulations show that in the next cou- 
ple of years the precision of these tests will increase by sev- 
eral orders of magnitude. The combination of several such 
systems in a PFE antenna array for the detection of PFE 
effects can be used to obtain a full sky coverage. 
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